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Multi-Dimensional Real Fourier Transform 


The problem: 

To compute the one-dimensional and multi-dimen- 
sional Fourier transforms for real data, the multi- 
dimensional complex Fourier transforms, and the 
multi-dimensional sine, cosine, and sine-cosine trans- 
forms. 


The solution: 

Four subroutines, each of which performs the 
required calculation for one of the four problems 
given. Each subroutine uses the Cooley-Tukey Fast 
Fourier Transform (FFT), and, in all but the one- 
dimensional case, the transforms are calculated in up 
to six dimensions. 


How it’s done: 


The Cooley-Tukey algorithm is an efficient com- 
putational procedure for evaluating expressions of 
the form: 
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where xj and a^ are complex, and N is composite. 
If the calculations were to be carried out in a 
straightforward manner, with an operation defined as 
one complex multiplication and addition, a total of 
N(N- 1) operations would be required; i.e., N-1 opera- 
tions for each value of j. The Cooley-Tukey al- 


gorithm, however, organizes the calculations in such 
a way that only (N log 2 N) operations are required, 
which represents a significant savings for large N. 

Notes: 

1. These programs are written in FORTRAN V for 
use on the UNI VAC-1 108 computer. 

2. The specific subroutines are as follows: NPO- 
11648, Multi-Dimensional Real Fourier Trans- 
form; NPO- 11649, One-Dimensional Real Fourier 
Transform; NPO-11651, Multi-Dimensional Com- 
plex Fourier Transform; and NPO-11652, Multi- 
Dimensional Sine, Cosine, and Sine-Cosine 
Transforms. 

3. Inquiries concerning these programs should specify 
the program desired, and may be directed to: 

COSMIC 
Barrow Hall 
University of Georgia 
Athens, Georgia 30601 
Reference: B71-10133 

Patent status: 

No patent action is contemplated by NASA. 

Source: Fred T. Krogh of 
Jet Propulsion Laboratory 
under contract to 
NASA Pasadena Office 
(NPO-11648, 11649, 11651, 11652) 
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